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ABSTRACT 

Recent analyses of datasets acquired at the Brookhaven National Laboratory and at the 
Physikalisch-Technische Bundesanstalt both show evidence of pronounced annual variations, sugges- 
tive of a solar influence. However, the phases of decay-rate maxima do not correspond precisely to the 
phase of minimum Sun-Earth distance, as might then be expected. We here examine the hypothesis 
that decay rates are influenced by an unknown solar radiation, but that the intensity of the radiation 
is influenced not only by the variation in Sun-Earth distance, but also by a possible North-South 
asymmetry in the solar emission mechanism. We flnd that this can lead to phases of decay-rate 
maxima in the range to 0.183 or 0.683 to 1 (September 6 to March 8) but that, according to this 
hypothesis, phases in the range 0.183 to 0.683 (March 8 to September 6) are "forbidden." We flnd 
that phases of the three datasets here analyzed fall in the allowed range. 

Subject headings: Methods: data analysis. Sun: particle emission 



1. INTRODUCTION 

iJenkins et all (|2009[ ) and iFischbach eFaLl (|2009D have 
drawn attention to the fact that decay-rate measure- 
ments made at the Brookh aven National Laboratory 
(BNL; lAlburger etUI (|1986l )) an d at the Phys ikalisch- 
Technische Bundesanstalt (PTB; iSiegert et ahl ([1998)) 
show clear evidence of periodic annual modulations. Ac- 
cording to their analyses, the decay rates appear to peak 
in winter, suggesting that some kind of radiation from the 
Sun is influencing the decay rates. However, the dates of 
maximum rate differ signiflcantly from January 3, which 
is the date of minimum Sun-Earth distance, which argues 
against a simple model in which decay rates are solely in- 
fluenced by an unidentifled but isotropic radiation from 
the Sun. The purpose of this article is to propose an 
interpretation of the phases of the annual oscillations. 

The suggestion that nuclear decay rates may be vari- 
able has of course been questioned by a number of other 
groups: 



1. iSemkow et al. I (|2009D and others suggested that 
these fluctuations have their origin in environmen- 
tal influences on the detec tor systems, or other 
systematic effects. However, IJenkins et al. I (|2010f l 
have shown that the results of their BNL and PTB 
analyses cannot be explained by variations of tem- 
perature, pressure, humidity, etc. 



2. iNorman et al. I ()2009t ). re-examined data from sev- 
eral nuclear-decay experiments, but found no ev- 
idence for a correlation with Sun-Earth distance. 
However, we have re-analyzed data for ^^Na and 
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^^Ti decay rates, which Norman and his collabora- 
tors generously provided to our collaboration. Nor- 
man and his colleagues compared these data with 
a sinusoidal modulation of fixed amplitude 0.15%, 
and fixed phase of year (0.01), to con form with the 
results found bv IJenkins et al.l (|2009l ) in their anal- 
yses of BNL and PTB data. However, we have 
carried out a power spectrum analysis of the Nor- 
man data, which searches for evidence of periodic- 
ity and yields esti mates of the amplitud e and phase 
of any oscillation to'Keefe et al. 1120111) . Our anal- 
ysis yields evidence of an annual modulation, with 
significance estimated at 1%. We estimate the am- 
plitude to be 0.034% [smaller than that considered 
by Norman et al. by a factor of 5], and the phase 
to be 0.06 ± 0.05, which is compatible with the 
value [0.01] assumed by Norman et al. The fact 
that the amplitude inferred from the data of Nor- 
man, et al. is smaller than what is observed in the 
BNL or PTB data is compatible wit h the standard 
weak interaction theory, as noted bv lJenkins et al.~1 
([2010). 



3. iCooped (|2009( ) analyzed data from the power out- 
put of the radioisotope thermoelectric generators 
(RTGs) aboard the Cassini spacecraft, but found 
no signifi cant deviations from exponential decay. 
However, IJenkins et al. I ()2010D have shown that 
there is in fact no conflict between Cooper's re- 
sults and their results. This is due in part to the 
fact that the Cassini RTGs derive their power from 
the alpha decay of ^^®Pu, whereas the periodic ef- 
fects observed in various data sets always involve 
beta decays. 

We have recently presented evidence that both BNL 
and PTB measurements also exhibit a periodicity sug- 
gestive of t he influence of solar r otation-possibly of the 
solar core (iSturrock et al.|[2"010al lH). Such a periodicity 
would indicate that the radiation is anisotropic, indicat- 
ing speciflcally a longitudinal (East- West) asymmetry. 
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If the radiation has a longitudinal asymmetry, it is rea- 
sonable to suspect that it may also have a latitudinal 
(North-South) asymmetry. 

The goal and scope of this article are limited to investi- 
gating the consequences of the following two hypotheses: 

(a) Some nuclear decay rates are influenced by an un- 
known solar radiation; and 

(b) The flux of this radiation, as it might be measured 
at Earth, is influenced by two factors: the variation 
in the Sun-Earth distance due to the eccentricity of 
the Earth's orbit, and a North-South asymmetry in 
the solar radiation. 

In Section O we analyze the properties of a model 
based on these two hypotheses. We find that, accord- 
ing to this model, the maximum decay rate should occur 
in the range approximately September 6 to March 8. In 
Section [31 we compare this model with data from the 
BNL and PTB experiments, and find that the phases of 
maximum variability appear to be compatible with this 
model. We discuss these results in Section H) 

2. MODEL 

It is convenient to measure phase of the year so as to 
run from to 1. We denote by 0o the phase of the min- 
imum Sun-Earth distance, and by 0a a phase related to 
the solar axis of rotation. Since the Sun-Earth distance 
is a minimum on January 3, 0o = 0.008. 

The Earth has maximum exposure to the southern so- 
lar hemisphere on March 6 and maximum exposure to 
the northern hemisphere on September 8, corresponding 
to phases 0.178 and 0.687, respectively. The fact that 
the difference is not exactly 0.5 is due to the eccentricity 
of the Earth's orbit. Approximating the elliptical orbit 
by a circular orbit for the purpose of representing the 
North-South asymmetry, we adopt (j)A = 0.183, which 
gives maximum exposure to the southern hemisphere at 
phase 0.183 (March 8) and maximum exposure to the 
northern hemisphere at phase 0.683 (September 6). 

We also point out that we are implicitly invoking a 
third hypothesis: 

(c) Whatever region of the Sun is responsible for a so- 
lar influence on decay rates has an axis of rotation 
that is indistinguishable from that determined by 
observation of the solar photosphere. 

The eccentricity of the Earth's orbit leads to a variation 
in the flux of the hypothetical radiation influencing decay 
rates, normalized to mean value unity, given by 

Fo = l + Ccos[2^(0-0o)], (1) 

where C — 0.0334. The variation of the flux due to a 
North-South asymmetry is given by 

Fa^I- AC cos [2Tr{^-(^ a)], (2) 
where the "asymmetry coefficient" A is chosen so that: 

(a) If |A| = 1, the amplitude of the variation due 
to North-South asymmetry is equal to the ampli- 
tude of the variation due to the eccentricity of the 
Earth's orbit, and 



(b) If ^ > 0, radiation from the northern hemisphere 
exceeds that from the southern hemisphere (and 
vice versa). Since C is small (and if A is not large), 
the total variation in the fiux (retaining the normal- 
ization to mean value unit) is given by 

F=l + C {cos [27r ((/) - <j>o)] - A cos [27r (0 - 0^)]} ,(3) 

If the average fiux is responsible for an increase in the 
decay rate by the small "coupling coefficient" F, then 
the annual modulation of the decay rate R, normalized 
to mean value unity, is given by 

R= l + rC{cos [27r(0-(/)o)] - Acos [27r (0 - (/)a)]}(4) 

[We should note, however, that it is conceivable that the 
unknown radiation may act on some isotopes to reduce 
the decay rate. In this case the coupling coefficient F 
would be negative.] 

We find from inspection of this formula that, since 
14' A — 0o| < 0.5, the maximum value of R is to be found 
only in the range 4>a — 0.5 to (f>A or, equivalently, in the 
two ranges to (j)A and (j)A + 0.5 to 1, i.e. to 0.183 
and 0.683 to 1. The range (j)A to 0a + 0.5, i.e. 0.183 to 
0.683, is "forbidden." [For the case that the solar radia- 
tion acts to suppress the decay rate, these ranges would 
be reversed.] One may understand this result by noting 
that if we sum two sine waves, with phases 0i and 4>2, the 
peak will be found in the range (pi to 02 if |0i — 02| < 0.5, 
and outside that range if |0i — 02| > 0.5. 

Expressed in the form in which it would be measured 
experimentally, this formula becomes 

i?= 1 + Xcos[27r(0-K)], (5) 

in which K is the amplitude of the variation, and k is 
the phase of the peak. By comparing Eqs. [5] and [S] and 
separating out coefficients of cos (27r0) and sin (27r0), we 
find that 

Kcos {2TTK)=rC [cos (27r0o) - A cos (27r0A)] 
i4:cos(27rK)=FC[sin(27r0o) - ^ sin (27r0A)] (6) 

from which we find that the asymmetry coefficient A is 
related to k by 

^ ^ sin [27r {k - 0^)] 
sin [27r {k — 0a )] 

This coefficient is shown as a function of phase in Fig- 
ure [TJ The coefficient is negative for k = to 0a and 
positive for k = 0a + 0.5 to 1. If k = 0o, the coefficient 
is zero, as we would expect. If k = 0a, the coefficient is 
infinite, again as we would expect. 

We also find from Eq. |6] that the coupling coefficient 
F is given by 

T = {K/C)G{t,), (8) 
where the "coupling factor" G (k) is given by 

^ ^ sin [2tt [k - 0a)] 

sin[27r(0A-0o)]' ^ ' 

This factor is shown as a function of phase in Figure 
[2l If K = 0o, the factor is unity. If k = 0a or 0a -I- 0.5, 
the factor is zero. 




Phases of Decay Rate Variations 

^ I 4>L,4>u denote the "l-sigma" offsets. 
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Fig. 1. — Asymmetry coefficient (A) as a function of phase ((/>), 
showing the values for BNL ^^Cl, BNL '-^^Si, and PTB ^26^^. 
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Fig. 2. — Coupling factor (G) as a function of phase (0), showing 
the values for BNL 3^01, BNL ^^Si, and PTB ^26^^. 

Figures [T] and [2] confirm that there are no valid so- 
lutions to the relevant equations for k in the range 
(t>A < n < 4'A + 0.5. 
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Fig. 3. — Power (5) as a function of phase {<j>) for BNL ^^Cl data, 
for frequency 1 yr~^. 
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Fig. 4. — Power (5) as a function of phase (0) for BNL •^2gi data, 
for frequency 1 yr^-"^. 



3. DATA ANALYSIS 

We now compare this model with BNL and PTB data. 
In our earlier power spectrum anal v sis ([Sturrock et alJ 
l2010aD . we followed lAlburger et al.l (|1986f ) and consid- 
ered only the ratio of the ^^Cl and ^^Si decay rates. 
However, since we are now interested primarily in the 
phases of the decay-rate oscillations (rather than their 
amplitudes), it is necessary to consider the two datasets 
separately. 

We determine the phase of the maximum in the annual 
variation of each dataset by using a likelihood method 
previ ously introduced for the analysis of solar neutrino 
data (jSturrock et al.l |2005() . However, rather than scan 
the power 5 as a function of frequency, we now set the 
frequency at 1 yr"-"^ and scan the power as a function 
of phase (f). We show the results in Figures [H IH and [5] 
for BNL 36C1, BNL 32Si, and PTB ^^^Ra., respectively. 
The results are summarized in Table [1] where we list, 
for each element of each experiment, the phase 0p of the 
maximum power, the maximum power Sp, the amplitude 
of the modulation, and the phases 0^, 0^/, for which the 
power is less than the maximum by 0.5. The phases and 
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Fig. 5. — Power (S) as a function of phase ((/■) for PTB '^'^^Ka, 
data, for frequency 1 yr~^. 

For '^^Si and ^^^Ra, we can state with some confidence 
that the phases fall in the allowed range. However, we 
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TABLE 1 

Phase, power and amplitude estimates of annual decay-rate variations in the BNL and PTE experiments. 



Experiment 


Eicmcnt 


0P 


6'p 


Amplitude 


<t>L 


'Pu 


BNL 




0.695 


14.46 


5.29x10"" 


0.688 


0.722 


BNL 


32 Si 


0.899 


3.52 


2.72x10-4 


0.836 


0.962 


PTB 


Ra 


0.085 


495.05 


8.40x10""' 


0.080 


0.090 



see from Figure [3] that, for ■^^Cl, we can state only that 
there appears to be 50% probabihty that the phase is in 
the allowed range. An interesting fact to note is that the 
■^^Si and ^^Cl were measured concurrently, in an alter- 
nating fashion, on the same instrument. The fact that 
these two nuclides do not show the same phase proper- 
ties further supports the consideration that the observed 
effects are not systematic in nature, as well as reinforc- 
ing the notion that different isotopes will have different 
inherent sensitivities. 

4. DISCUSSION 

It is important to note that our model refers implicitly 
to the timing of the unknown solar radiation responsible 
for the decay-rate variations. Since there may be a time 
lag between the arrival of the radiation and the emission 
of a beta particle, there may be a discrepancy between 
the timings of the two events. This point is relevant to 
the BNL measurements of the '^^Cl and '^^Si count rates 
summarized in Table [1] 

What is experimentally observed in "^^Si decay is the 
1709 keV beta particle emitted by "^^P in the sequence 
32gj _^32 p ^32 gjjjj,g ^j^g 225 keV beta particle emit- 
ted from the transition has too low an energy to have 
been detected by the BNL experiment. Since the half- 
life of ^^P is 14.28 days, '^^P comes into equilibrium with 
the parent "^^Si on a time scale of order 14.28d/ln (2). 
Hence the response of the beta-detector lags the initial 
^^Si decay by about 21 days. By contrast, ^^Cl decays 
directly to ^^Ar with the emission of a 709 keV beta 
particle. Hence the calculated phase difference between 
the '^^Si and ^^Cl measurements of about 75 days, shown 
in Table [U is probably due in part to the mechanisms 
responsible for the decay-rate variations. 

More generally, we draw attention to the following 
points: The phase distribution that follows from our 
model rests on the assumptions that decay-rate varia- 
tions are due entirely to an unknown solar radiation, and 
that the axis of rotation of the core is parallel to that of 
the photosphere. If we find that, for some isotopes, the 
timing of decays does not conform to the expectations 
of this model, we must question these assumptions. One 
would then need to consider the possibilities that decays 
may be due in part to a cosmolo gical particle o r radia tion 
flux, as has been suggested bv lBaurov et al.l ()2009 f). or 
that the axis of rotation of the solar core is not parallel 
to that of the photosphere. 

As we see from Table [1] there is a huge variation in 
maximum power among the three cases, but a compar- 
atively small range in amplitude. The very large power 
for the PTB case is due in part to the longer duration 
of the experiment (14.8 years for PTB as against 7.8 
years for BNL). However, even allowing for this differ- 
ence, we see that different elements have very different 
sensitivities to whatever is causing the annual modula- 
tion. As noted above, this is consistent with expecta- 



tions from standard w eak interaction theory, as noted by 
iJenkins et alTI (|20Tot) . 

For the purposes of this article, we refer to "decays" 
as if a decay event is a simple process, but this is typi- 
cally not the case. For instance, the decay of ^^^Ra is a 
complex process. Although the decay of a ^^^Ra nucleus 
begins with alpha emission, and although after about 20 
years ^^^Ra nuclei are in secular equilibrium with the 
majority of its 13 radioactive daughters, it takes about 
200 y ears for the pro cesses to effectively reach equilib- 
rium (IChristmas et al.1119831 : IChiste et al.1[200l . When 
equilibrium is reached, 42% of the photon emissions are 
due to beta-decaying daughters. The ionization chamber 
utilized in the PTB experiment cannot discriminate be- 
tween either (alpha or beta) type of decay: the chamber 
measures only the total energy deposited by the incident 
photons, which have their origins in both types of decay 
from several different isotopes. 

The asymmetry coefficient and the coupling factor, 
given by Eqs. [7] and [9l respectively, are listed in Ta- 
ble [T] and shown in Figures [T] and [21 respectively. We see 
that the phases fall in the "allowed" range 0.683 to 0.183 
(avoiding the "forbidden" range 0.183 to 0.683). We 
identify in these figures the phases of maximum power 
listed in Table [T] On applying the results shown in these 
two figures to the data shown in Table [U we obtain the 
values of the asymmetry coefficient and of the coupling 
coefficient listed in Table [2] 

In reviewing the results shown in Table [2l we see that 
the results for the ^^Cl experiment are ill-determined, 
since the phase is close to the critical phase 0.683. In 
comparing the results for the BNL ^^Si experiment and 
the PTB ^^Ra experiment, we see that there is no sub- 
stantial difference between the estimates of the magni- 
tudes of the asymmetry and coupling coefficients, de- 
spite the fact that the powers differ greatly. It is inter- 
esting to note the difference in sign between these two 
experiments: for both elements investigated in the BNL 
data, the modulation shows evidence of a larger influence 
arising from the northern solar hemisphere, whereas for 
the PTB experiment, the modulation shows evidence of 
a larger influence arising from the southern hemisphere. 
Since the time intervals of the two experiments are differ- 
ent (1982.13 to 1989.93 for BNL and 1983.86 to 1999.17 
for PTB), we cannot infer whether the difference is due 
to variability of the solar influence, or to a difference in 
the "sensitivity" of the specific decaying isotopes. It will 
therefore be important to study a variety of isotopes con- 
currently for precisely the same time interval. A frame- 
work to descr ibe a possib le class of mechanisms is pre- 
sented in Fischbach et al.l (12009). 

We h ave found evidence, r eported in previous publi- 
cations (jSturrock et al.|[2"010a ibV of modulations in de- 
cay rates that may be attributable to rotation at or 
near the solar core. We have also recently found ev- 
idence of a Rieger-type modulation that we may at- 
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TABLE 2 

Estimates of the asymmetry coefficients (A) and coupling coefficients (F) derived 

FROM DATA FOR THE BNL AND PTE EXPERIMENTS. 





TV,1 OTTl Oirf" 


L/il Lr LI. Lit; 




Asymmetry 

V^tJtJliiLlL-iiL 


Oouplmg 

V-'ULliiL-ll^liL 




36 CI 




0.688 


N/S 


N/S 


BNL 


5.29x10-"' 


0.695 


12.25 


0.0013 








0.722 


4.02 


0.0043 








0.836 


1.08 


0.0076 


BNL 




2.72x10-"' 


0.899 


0.65 


0.0090 








0.962 


0.29 


0.0090 








0.080 


-0.72 


0.0171 


PTB 




8.40x10-"' 


0.085 


-0.80 


0.0163 








0.090 


-0.89 


0.0156 



tribute to an "inner t achocline" separating the core 
from the radiative zone (jSturrock. Fischbach fc Jenkins I 
1201 If ). These results show that whatever radiation is 
responsible for influencing decay rates propagates essen- 
tially freely through the solar radiative and convection 
zones, suggesting that some flavor of neutrinos may be 
responsible for decay-rate variations. However, the stan- 
dard theory of neutrino physics does not provide for a 
mechanism by which neutrinos could have such an effect, 
raising the possibility of new physics. 

In view of these conflicting considerations, it is clearly 
desirable that we obtain further information relevant to 
the mechanism by which the Sun influences decay rates. 
If neutrinos or some other particles (or quanta of ra- 
diation) influence radioactive atoms, it is conceivable 



that they may transfer their momentum to the decay- 
ing nuclei, and thus exert a force on a bulk sample of 
radioactive material. This raises the possibility of car- 
rying out a "fifth force" type of experiment on a macro- 
scopic radioactive sample, using techniques recently de- 
veloped to search for deviation s from Newtonian gravity 
(jFischbach fc Talmadge I[l999l ). 

We are indebted to D. Alburger and G. Harbottle for 
supplying us with the BNL raw data, and to H. Schrader 
for supplying us with the PTB raw data. We are also in- 
debted to an anonymous referee for suggestions that ma- 
terially improved the article. The work of PAS was sup- 
ported in part by the NSF through Grant AST-06072572, 
and that of EF was supported in part by U.S. DOE con- 
tract No. DE-AC02-76ER071428. 
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